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ABSTRACT 


The  effect  of  the  wall  on  the  motion  of  elongated  aerosol 
particles  was  calculated  on  the  basis  of  the(powerful)  slender 
body  theory  for  a  quiet  and  flowing  medium.  The  cases  treated 
were  those  of  cylindrical  particles  moving  in  a  quiet  environ¬ 
ment  or  a  Couette  or  a  PolseUille  flow  near  plannar  surfaces  with 
various  inclinations  to  the  horizontal. 

Experiments  of  trajectory  photographing  were  performed  in  a  liquid- 
tank  under  creeping  flow  conditions  and  typical  (aerosol)  Reynolds 
numbers  whereby  very  good  correspondence  with  computations  was 
obtained. 

The  Brownian  diffusion  of  (very  small)  cylindrical  or  discoidal 
aerosol  particles  was  also  studied  while  being  based  on  Brenners' 
formal  structure.  To  evaluate  the  signficant  orientation  density 
function,  the  (mid-diameter) rotational  diffusion  coefficient  of  the 
particles  was  calculated  by  a  statistical-mechanical  method.  The 
values  of  the  coefficients  were  found  to  be  larger  than  those  computed 
on  the  basis  of  continuum  theory  by  a  factor  of  about  30. 

The  typical  diffusion  situation  of  deposition  on  a  plane  within  a  quiet 


air  was  treated. 
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INTRODUCTION 


During  the  past  several  years  a  series  of  studies  on  the  transport 
mechanics  of  non-spheric^.  aerosol  particles  was  carried  out  mostly 
with  the  support  of  the  European  Research  Office,  U.S.  Army. These  studies 
were  summarized  in  the  final  reports: 

The  Dynamics  of  Nonspherical  Particles 

I.  Controlled  Generation  of  Nonspherical  Aerosols  and  Methods 
for  Their  Size,  Concentration  and  Electric  Size  Measurement  (Contract 
DAJA  37-72-3875;  1974). 

II.  Mobility  in  Non-Orienting  Fields  (Contract  DAJA  37-74-C- 
1208;  1975). 

III.  Mobility  and  Deposition  in  Still  Air  (Contract  DAER0-75-G-021 ; 
1977). 

IV.  Deposition  from  Flow  (Contract  DAER0-76-G-014;  1978) 
and  in  the  parallel  journal  publications: 

1.  "On  the  stochastic  Motion  of  Nonspherical  Aerosol  Particles. I .  The 
Aerodynamic  Radius  Concept"  by  Isaiah  Galllly  and  Aaron  Hi-Cohen,  J.  Colloid 
Interface  Scl.  56,  443-59  (1976) 

2.  " _ ,11.  The  Overall  Drift  Angle  In  Sedimentation",  by 

Matitia  A.  Weiss,  Aaron-Hi  Cohen  and  Isaiah  Gallfly,  J.  Aerosol  Scl.  9» 
527-41  (1978). 

3.  "On  the  Orderly  Nature  of  the  Motion  of  Nonspherical  Aerosol 
Particles.  I.  Deposition  from  a  Laminar  Flow."  by  Isaiah  Galllly  and 
Alferd  D.  Eisner,  J.  Colloid  interface  Sci.68,  326-37  (1979). 
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Likewise,  they  were  reported  in  the  following  articles: 

"On  the  Orderly  Nature  of  the  Motion  of  Nonspherical  Aerosol 
Particles.  II.  Inertial  Collision  between  a  Spherical  Large  Drop¬ 
let  and  an  Axially  Symmetrica)  Elongated  Particle"  by  Isaiah  Gal  lily 
and  Aaron-Hi  Cohe,  J.  Colloid  interface  Sci.  68^  338-56  (1979). 

5.  " _ ,11. _ ,  Further  Results,"  by  Isaiah  Gallily  and 

Aaron  Hi-Cohen,  ibid.  7J_,  628-30  ( 1 979)  • 

The  efforts  invested  in  ail  of  these  investigations  culminated  in  the 
establishment  of  a  method  of  claculation  of  the  gravitational  sedimentation 
in  still  and  laminarily  flowing  gas,  and  of  the  inertial  impaction  on 
various  surfaces  of  ax i -symmetrical  particles  such  as  cylinders  and 
disks.  The  method  has  been  el laborated  in  the  (above)  publications  1-5  and 
partially  checked  by  experiments;  however,  it  is  far  from  being  complete. 

One  problem  in  the  analysis  of  our  particles'  transfer  phenomena  which 
has  not  been  considered  is  the  effect  of  the  surface  exposed  to  the  aerosol  on 
the  collection  process  itself.  The  effect,  related  here  only  to  the  motion 
of  the  particles  in  the  adjoining  gas  layer,  is  due  to  ; 

a.  The  fluiddynamic  interaction  between  a  particle  approaching  a 

surface  and  the  latter  which  essentially  stems  from  the  change  of  the  flow  field 
due  to  the  emergence  of  an  additional  boundary  condition. 

b.  The  presence  of  long  (electrostatic)  or  short  (Van  der  Waals)  surface 
particle  forces. 

Another  untreated  problem  concerns  the  stochastic  dfffusional  turbulent  and/or 
Brownian  movement  of  the  particles. 
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It  is  needless  to  stress  the  significance  of  the  first  problem  which  has 
been  treated  scarcely  for  the  case  of  spherical  particles  by  Brenner(l), 
Gal  lily  and  Mahrer  (2)  and  others  and  has  been  almost  untouched  In  the 
case  of  the  nonepherical  ones. 

In  most  studies  on  the  orderly  motion  of  any  particle,  expressed  by  the 
equation  of  translation 


ynf/tt/c/-f=^Fe:+Fr  W 

t 

and  rotation 

d(I*.c*)/dt  =Z  Le+  Lf  [2] 

the  fluiddynamic  operating  force  Fp  and  torque  Lp  have  been  calculated 
for  the  substitute  situation  in  which  the  particle  is  infinitely  far  from 
the  (col lect ing) surface .  This  has  been  obviously  a  simplification  whose 
justification  was  checked  in  the  calculation  of  Fp  and  Lp  for  some  spheri¬ 
cal  cases  (1,2).  For  the  nonspherlcal  case,  it  has  not  been  studied  at  all. 

The  importance  of  the  second  problem  is  self  evident  for  sufficiently 
turbulent  field  or  very  small  particles  for  which  the  fluctuating  part  of  the 
force  Fp,  viz.  Fp  -  <Fr)  »  can  not  be  neglected  In  comparison  with 
the  mean  one  <Fp>. 
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AIM  OF  STUDY 

In  view  of  the  above  review  of  the  state  of  research,  it  was  thought 
worthwhile  to  investigate: 

1.  The  effect  of  the  collecting  surface  on  the  motion  of  nonsphe- 
rical  particles  near,  and  their  deposition  on  it,  and,  independently, 

2,  The  diffusion  of  nonspherical  particles. 

These  two  problems  seemed  immediately  too  formidable  for  a  general 
treatment;  thus  their  study  was  restricted  in  the  first  subject  to  the 
case  of  cylindrical  particles  much  larger  than  the  mean  free  path  1  of 
the  medium's  molecules  (  K-vi  j  ),  and  a  planar  solid  wall,  and 
in  the  second  subject  to  the  Brownian  diffusion  of  cylindrical  and  discoldal 
particles  much  smaller  than  l  (  <<  *  )  ,  Cylindrical  and 

discoidal  forms  respectively  serve  as  idealized  model  shapes  of  the 
real  (harmful ) asbestos  or  other  significant  fibrous  particles  and  of 
"1atelet-1 ike  ones.  So,  the  restriction  to  the  above  morphology  is  not 
very  limiting  in  significance. 


A.  THE  WALL  EFFECT 


I.  Theoretical 
1 .  Method 

a.  General 

Here  we  preferred  to  use  slender  body  theory  rather  than  the  (limited) 
method  of  Wakiya  which  is  based  on  Fourier  analysis  technique  (3). 


The  slender  body  theory  method  is  much  more  general  than  the  latter; 
it  can  be  employed  in  situations  where  the  (aerosol)  particle  is  even 
far^rom  walls  with  a  better  advantage  than  that  obtained  by  being  based 
on  the  much  used  Oberbeck's  (^)  and  Jeffery's  (5)  equations.  In  other 
words,  it  is  very  useful. 

The  method  is  apt  for  the  treatment  of  the  motion  of  particles  whose 

characteristic  size  r  c  (=  V  &  /  M  n  )  is  much  larger 

—  ,  .  7 

than  the  mean  free  path  1  of  the  mediums  molecules  (  K'-y,  -  —  »  3 

re 

and  whose  characteristic  Reynolds  number  Re  is  much  smaller  than 
unity.  In  fhat  case,  the  equations  of  fluid  motion  to  be  solved  are 


[a  V  vf**)  -  Vf  r  0 


[31 


and 


7 •  v  (r)  = 


M 


The  theory,  developed  by  Batchelor  (6)  ,Chwang  and  Wu(7) ,  Blake  (8), 

Cox  (9)  and  others,  deals  with  elongated  particles  in  which  the  ratio  fc 
between  the  characteristic  cross-sectional  radius  b  and  the  semi-length  1 
is  W  =  -S-  «  i  ,the  change  of  the  local  equivalent  cross-sectional 

'For  symbols  see  NOMENCLATURE 


radius 


kr,  is  small  compared  with  (rr<  /  (r 
of  curvature  is  much  greater  than  t  (Fig.  1). 


and  the  local  radius 


P' 


Fig.  1.  A  slender  body  (not  necessarily  a  straight  cylinder); 
coordinates  used  in  the  analysis. 


Under  these  coidi tions ,  the  disturbance  velocity  v  (r)  produced  by  the 
moving  particle  is  given  by  the  multipole  expansion 
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e  ( 

v  1  C*" )  =  (  y  O')  u  ij  (.>,  r'}  otr*  +  j"  ijk  (*■')  Uijk  Cr,  r‘>  3I  r*1 

-t  -t 

+  T-fj).*.  O'")  O,  *•')  ofr'  + 


where,  with  the  Einstein's  convention  adhered  to,  (r>)  (Stokeslet), 


•fjk  c  *-' ) 


(doublet)  , 


■f  i  10  C  **' ) 


etc.  are 


(imaginary)  poles  causing  the  disturbances  ay  (r,  r‘)  etc,;  the  latter 
are  actually  the  proper  Green  functions  (APPENDIX  1). 


The  calculation  of  the  (critically  required  )  fluiddynamic  force 
and  torque  U  f  appearing  in  Eqs.  [1]  and  [2]  can  be  performed 
once  the  Stokelet  -f.  (r1)  distribution  is  known  since  (Flg.l) 


ff  =  [  f  CO  cli 


uF  =  J  (  x  x  -f  (*)  )  c(s  1 7 ] 

-l 

The  finding  of  (**')  cannot  be  carried  out  by  the  solution  of  the 

full  Eq.  [5]  with  the  proper  boundary  conditions  on  the  particle 
(v  (r)-  V  (r)  **t  U  (r)  +  w  x  x  for  *■  — *  *  +  b)  and 

on  the  wall  (v(tr')  =  0  )  •  Thus,  one  uses  as  an  approximation  for  t  »  fr 

(Fig . 1 )  an  "outer  expansion"  in  which 


I 


8 


l 

v;  O')  Ujj  (r,r')  dr1  C83 

-l 

only,  and  for  4  <*'.  ;  <*■  l  (for  which  the  disturbance  can  be  assumed  to 
evolve  from  a  portion  of  an  infinite  straight  cylinder  of  radius  4r, 
an  "inner  expansion"  (6) 


<«o 


_j _ 

2rr  u 


O)  [U 


+  K(«)]  +  Uv(») 


C91 


and 


vf  Cr)  =  -  C  Cr)[6n  —1*  +  -+  l<u  (*)J  +  U-t» 

4  TT  u.  ?  o1  v 


where  l1  ,  L'  (  =  a  ,  3  )  respectively  stand  for  the  particle-locked 

local  coordinates  x,.  ,  xx.  ,  x,f  ,  C,',  ,  C  •  are  quantities 

affected  by  the  wall  and  K  ,  Kt  •  are  constants  defined  by  the 

c 

geometry  of  the  cross-section  at  e'  (Fig.  1).  This  "Inner  expansion" 
can  be  also  written  In  the  vectorial  form  fit  for  any  coordinate  system 


4irpvO)  =  M  •  0  bi  -y*  +  M  •  TT  •  C  +  I*C-I*C  +  ^tt/aU  [tl] 

where  H  ,  C,  TT,  I,  I'  are  defined  In  APPENDIX  1. 

The  two  approximations  are  matched  in  a  region  in  which  they  are  both 
assumed  to  hold  (6-«  ?  )  .  and  the  relevant  constants  so  derived. 


In  the  case  of  the  "outer  expansion",  one  can  write 


9 


u  ;  (w^r') 


O,*' 


(  r,  •") 


Cl  2] 


| 

where  u;j  (r,  r'j  has  singularities  on  the  particles  axis 

W 

and  U-y  (  rf  r' ,  is  non-singular  when  §  »  £ 


b.  Straight  Body 

In  this  case  Eq.  [8]  is  integrated  by  Batchelor's  approximation  to  [11}  ((,) 

which  gives 

kitfL  v(r)  ={-£-  +  - 

+  —  M  •  R  +  A_  +  I*-f  [13] 

I  "  ,  R  and  .A.  being  defined  In  APPENDIX  1  ,  <T  =  [  tr »  (2  £/(>■)] 

and  X  *  (^r1  /  &)  .  Now,  matching  the  coefficients  of  the 

(singular)  fr?  part  in  [11]  and  [13]  one  gets 

t  * 

c(«)  =  f(«)  ,  [14] 


and,  doing  the  same  for  the  non-singular  part, 

tticp  U(»)  - ('~gr  ♦  f*i[  ^  W*f  +  ^  'H'4+-jpM*R+A  PS) 

where  1  °  =  -l"  +  \  l'  and  Tt  *  =  TT  •  H  . 
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At  this  stage  Cl)  can  be  jjjB.rlc.llY  calcelated  free,  [151 
by  dividing  oar  particles’  axis  Into  H  egaal  parts,  to  each  of  which  Is 

associated  a . .  f ,  Cs)  .  «*  b-for»ln,  the  miration  of 

(  r  (  r1  )  and  — respectively  appearing  in 
JW  !nd  B  •,  so,  for  the  two-dimensional  case  of  particles'  motion, 
one  obtains  2  N  equations  for  2  H  unknowns  while  for  the  three- 
dimensional  case  one  obtains  3  *  equations  for  3  »  unknowns. 

According  to  Batchelor  W,  It  is  possible  to  use  here  a  quick  iterative 
method  bv  writing 


D*3 


1  .The  Wall  Effect 


Dealing  with  a  cylindrical  particle  undergoing  a 
and  rotational  motion  near  a  planar  solid  wall  (Fig. 


translational 
2)  and  statisfylng 


n 


Fig. 2.  A  slender  cylinder-planar  solid  wall  system,  coordinates 
used  In  analysis  , 


the  no-sllp  condition  on  all  boundaries,  one  has  to  Introduce  on  the 
reflected  Image  of  the  particles'  axis  a  l Ine-dlstrlbutlon  of  point 
forces,  -fa(x)  ,  a  1 1ne  distribution  of  doublets,  Jk-fjOr)  • 
and  a  line  distribution  of  quodrupdes,  -  (*)  ^  • 
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In  this  case 


£  Li  - 

Ir-VM 


,C  r,  -  r,’)  ( 

Ir-rM* 


[173 


and 


U,  ,N  JV_  .  Jr<  -  ^  {yJl  rA-+*^  <*)  (fy* *<* k-  -  h  5  ' 
^r»  r  '  “  lr'1  r  ; 


vs.u  aLn%;  f>*  r,%  [,8] 

"  lr*M 5  * r  r 

,  ,  afMl  l  i  k  -  1,2*3  respectively  . 

where  *  i»  1  or  2  and  SJ  , 


i  1 ,  Mot  1  on_w jth t rj_a_F J,ow 

•**"  expansion"  In  a  onl.t  >«  ’»  ’~,H* 

to  write  In  this  case  for 


v’(r)-[\  Cr)Uy(r,r->otr' 

-l 


P91 


v-W  .s,H. 

undisturbed  flew  ((v<r)  V*(r)  **  r  *  ~ 

expansion"  can  b.  express.*  W  »  *°  M 

*  (  r  «  C%  (V  )  I*  now  affected  by  the 

but  where  G  55  \  ''  i  ?  *  *  * 


flow. 
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Thus,  after  performing  the  necessary  matching  between  inland  the 
modified  [It.  one  gets  like  PS 


-  TT>  N>-f  s' 


M  .  R  +  .A.  (*; 


[20] 


in  which 


U1  («)  =  Cl  (si  -  v‘ 


< 


U(f'  is  defined  by 


v(r)  — f  «ri  for  r 
as  marked  above. 


+  b  , and  ail  other  symbols  are 


1  axis  in  the 


Expanding  V 
fora 


*(r)  on  the  (straight)  particles’  ax 


M 


and  putting  x K=  x ^  (i4  being 
one  has  for  U  Cs) 


oos  6  and  i,  being  sin^  )> 


u 


!  (x)  =  Ut(x) -vT(x)  a  Bi,X  -3^ 


pi] 


.  Is  the  rotational  velocity 
where  B;,  —  £  W  lw  ”  *•  }  L 

of  the  particle  aod  V;  Is  Its  translation!  velocity. 


t 


V 
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2.  Computational  Procedure 

The  trajectlons  of  our  particle  were  derived  from  the  solution  of  Eqs.  [1. 
and  ~JZ]  on  the  specific  assumptions  that  the  particles  are: 

1.  homogeneously  dense 

and 

ii.  affected  only  by  the  gravitational  force  mj,  the  fluiddynamic 
force  Fp  (Eq.  [6]  )  aid  the  fluiddynamic  torque  Lp. 

At  first  were  treated  two-dimensional  cases,  in  which  particles  axis 
stays  in  a  plane  containing  g  ;  then  we  considered  three  dimensional 
situations. 

For  the  two-dimensional  cases,  where  the  [I]  and  [2]  turn  into 


wc(V//c(i  =  (-W1--WI  ')3  +  F  F 

and  (since  cj,  -  coa  =  o) 


I  -j  ^5  —  ^3  * 


[23] 


[2*] 


we  employed: 

i.  A  quick  technique  based  on  Eq.  [16] 
and 

11  An  elaborate  technique  described  In  APPENDIX  1. 

The  computation  itself  was  performed  (on  the  Hebrew  University  CYBER  7^  CDC 
puter)  by  the  use  of  a  fourth-order  Range  -Kutta  method;  the  time  intervals 
were  from  5  •  IQ-1*  to  3  *  10-1*  second,  depending  on  the  (decreasing) 


corn- 


gap  between  the  particle  and  the  wall. 
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II.  Experimental 

The  experiments  were  essentially  intended  to  check  the  calculations  in 
their  significant  points.  Thus,  one  could  in  principle  either: 

1.  Determine  the  force  Fp  and  torque  Lp  from  measured  trajectories  and 
orientations  and  compare  with  theory 

or 

2.  Determine  the  whole  trajectories  of  the  particles  with  their  orienta¬ 
tions  and  test  against  computations. 

At  any  case  it  was  necessary  to  measure  particles'  location  and 
orientation  as  a  function  of  time;  so,  since  the  observation  of  small  cylind¬ 
rical  aerosol  particles  is  complicated  by  light  scattering  phenomena 
(10),  it  was  thought  worthwhile  to  use  large  (macro)  particles  which  move  in  oil 
in  a  creeping  motion  (  Re  <<  ±  )  too. 

1 .  Method 

In  essence  and  without  loss  of  generality,  our  experiments  were  carried 
photographing 

out  by/&  two-dimensional  motion  of  a  homogeneously  dense  cylindrical  body 
falling  under  the  action  of  gravitation  in  an  oil  medium.  Here,  out  of  the 
above  two  alternatives,  we  chose  the  latter  which  is  more  exact  than  the 
former  as  it  does  not  involve  differentiation  of  location  and  velocity  data 
curves. 

2.  Apparatus 

The  apparatus  designed  and  constructed  for  the  above  purpose  (Fig.  3, a 


and  b)  consisted  of: 


Fig.  3.  The  experimental  apparatus  (achematlcal )  (a)  A,  Tank*, 

B,  plate  (wa!1);C,  release  device;  D,  cylinder;  E,  movie-camera;  F,  platform; 
G,  scale,  H,  mirror,  (b)  K,  cylinders'  holding  part  of  release  device; 

L,M,  sliding  rod  and  guide,  respectively,  a^,  82  pivote  rods. 
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A  tank  A  containing  silicone  oil  in  which  was  suspended  a  polished  steel 
plate  B,  a  device  C  to  release  a  metallic  cylinder  C  Into  the  liquid,  a 
movie-camera  E  and  light  source  for  the  trajectory  photographing,  a  moving 
platform  F  for  the  camera,  and  a  scale  G  -mirror  -  light  source  combination 
for  locating  the  particle. 


The  tank  was  made  of  Plexiglas  in  the  (inner)  dimensions  of  50  cm  height, 

50  cm.  breadth  and  3  cm.  width  and  was  kept  vertical  by  Its  own  weight. 

The  liquid  (F  111/10000  dimethyl  silicone  fluid,  ICI,  England)  had  a 
(nominal)viscosity  ^  of  10^  centistokes  or  about  100  poise  dynamic 
viscosity  at  25°C  with  a  temperature  coefficient  of  0.02  %  /°C,  and  a 
density  of  0.968  g.  cm  ^  at  25°C.  This  11  quid  was  satisfactorily  Newtonian 
as  tested  before  each  experiments  by  measuring  p  at  various  shear  rates 
3  V;  /  9  by  the  settling  velocity  of  various  sized  Stokesian  steel 

balls.  The  cylinder  used  as  a  tested  particle  had  a  radius  b  of  0.0075  cm, 
a  semi-length  1  of  1.5  cm.  and  a  density  of  7-780  g.  cm  It  was  intro¬ 
duced  into  the  liquid  by  a  slow  sliding  movement  in  the  groove  of  the  devices' 
edge  K  (Fig.  3  b)  and  along  the  middle  (vertical)  line  of  the  plate 
(wall).  The  camera  was  a  Bolex  H  16  Reflex  movie  instrument  (Bail  lard, 
Switzerland)  equipped  with  a  1:  18  lens. 
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3.  Procedure 

The  viscosity  p.  of  the  liquid  was  determined  before  every  run  as 
described  above;  it  was  found  to  be  constant,  within  the  experimental 
error,  at  the  value  of  112+2  poise  for  16. 4°C .  The  temperature  of 
the  oil,  measured  at  various  locations  in  the  tank,  was  constant  to  a 
degree  of  0.02°C,  which  amounted  to  a  viscosity  change  of  0.02  poise. 

The  trajectory  photographing  was  performed  while  keeping  the  center  of 
the  falling  cylinder  on  the  optical  axis  of  the  camera  as  well  as 
possible.  The  reading  of  position  and  orientation  of  the  particles  image  on 
the  movie  film  was  carried  out  with  a  Nikon  6CT2  profile  projector  (Nikon, 
Japan)  to  a  (nominal)  accuracy  of  +  0.001  mm. 

Here  it  should  be  pointed  out  about  the  critical  requirement  of  two- 
dimensionality  to  which  a  self-restriction  was  imposed.  Due  to  the  cylindri¬ 
cal  rod  release-operation, there  was  always  a  chance  that  the  particle  will 
not  move  with  its  axis  in  the  plane  of  g  .  So,  in  every  experiment,  the 
film  image  of  the  falling  rod  was  aposteriori  checked  against  the  expected 
value  deduced  from  the  actual  size  and  orientation  of  the  body,  and, when 
a  deviation  occured,this  run  was  discarded. 

III.  Results 

As  the  checking  experiments  had  to  be  restricted  to  a  liquid  medium 
and  macro  particles,  the  calculations  were  confined  to  that  situation  too. 
Oealing  with  physical  cases  which  can  be  analyzed  by  continuous,  classical 
fluiddynamic  methods,  this  did  not  matter  as  long  as  creeping  motion 
conditions  (  Re  <<  ±  x-n  <<  i  )  still  prevailed  and  the  liquid 

specific  effects  associated  with  the  acceleration  of  the  particles,  viz.  the 
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added  mass  and  moments  of  inertia  and  the  Basset  term  influence  (11), 
could  be  neglected. 

The  cases  selected  for  computation  were  chosen  to  be  a  two-dimensional 
motion  in  a  still  environment  near: 

i.  A  vertical  wall  (Figs.  5  through  8), 

ii.  An  horizontal  wall  (Figs.  9  through  11), 

iii.  An  inclined  wall  (Figs.  12  through  19) 

and  a  two-dimensional  motion  within: 

i.  A  two-dimensional  Couette  flow  over  an  horizontal  wall  (Fig.  20) 

ii.  A  two-d imens iona  1  PoiseuilJft  flow  through  an  horizontal  channel 
(Fig.  21). 

All  of  these  cases  simulate  real  situations.  Thus,  for  example,  if  the 
particles  are  small  enough,  then  any  corrugated  surface  can  be  represented 
by  inter-connected  plan  ar  portions  of  various  inclinations  to  the  horizontal 
likewise,  a  general  laminar  field  can  be  approximated  by  a  series  of  inter- 
adjoining  Couette  Hows.  The  restriction  to  two-dimensional  motion  which 
is  relaxed  latter  on  in  a  further  study  (not  reported  here)  is  not  serious 
since  in  many  instances  such  a  situation  is  indicative  enough  of  the  physical 
occurences . 

The  distances  and  angles  employed  in  the  calculations  are  brought  out 
in  Fig.  k  (the  unit  of  distance  is  always  the  semi-length  1).  Thus  the 
undisturbed  Couette  flow  was  expressed  by 


v/  Cx.)  =  k,  D, 


[25] 
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where  *,  is  the  (constant)  velocity  gradient,  while  the  (undisturbed) 
Poiseuille  flow  was  given  by 

(x.  ,  =  (  a2  -  D,s,  [26] 

where  is  a  (flow)  constant  and  a  is  the  half  width  of  the  channel. 

The  typical  results  shown  here  are  presented  in  terms  of  the  (operationa¬ 
lly)  important  drift  displacements  and  linear  velocities  and  the 
intuition-assisting  rotational  velocities  and  particle  orientation. 


Fig.  i*.  A  slender  body  near  a  wall.  Distances 
and  angles  used  in  calculations. 
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Fig.  7.  Orientation  and  angular  velocity  of  a  cylinder  near 
a  vertical  wall  as  a  function  of  time,  for  cases 
A, 8  of  Fig.  6. 
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Fig.  8.  Vertical  ( V  ^ )  and  hor  izontal^)  velocities  (in  l/sec) 
of  a  cylinder  near  a  vertical  wall  for  cases  A  and  B  of 
Fig.  6. 


Fig.  9.  Vertical  velocity  of  a  cylinder  falling  towards 
an  horizontal  plane  as  a  function  of  Dj  and  for 
various  initial  orientations  ©0  (in  rads). 

Same  cylinder's  and  liquids  properties  as  of  Fig. 5 


g.  10.  Drift  (horizontal)  velocity  of  a  cylinder  falling  towards 
an  horizontal  plane  for  various  Initial  orientations  Q0 
(In  rads).  Same  cylinder's  and  liquid's  properties  as  of  F 
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Fig.  11.  Angular  velocity  of  a  cylinder  falling  towards  an  horizontal 
plane  for  various  initial  orientations  (In  rods)  at  JD,  c 

*  10  1,  5  1.  Same  cylinders'  and  liquids  properties  as  of  Fig. 5. 


Fig.  12.  Time  history  of  a  cylinder  falling  towards  a  plane  inclined  to 
the  horizontal  at  angles:  A-15°,  B-30°,  C*A5°,  0-60°,E-75°; 
initial  orientation  -  vertical;  D,/0  •  10  1 .  Numbers  near 
positions  indicate  elapsed  time  in  seconds.  Same  cylinders  and 
mediums  properties  as  of  Fig.  5- 
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Fig. 


13.  Time  history  of  a  cylinder  falling  towards  an  inclined  plane. 
Same  parameters  as  of  Fig.  12  except  for  horizontal  initial 
or i entat ion . 


Fig.  19.  Angular  velocity  of  a  cylinder  falling  towards  an  inclined 
plane.  Same  parameters  and  conditions  as  of  Fig.  13.  ( 
reaches  a  minimum  at  the  point  P.). 


Motion  within  a  Flow 


Fig.  20.  Change  with  time  of  the  vertical  velocity,  A  (in  1/sec) .angular 
velocity  B  (in  rads/sec) ,  distance  from  the  wall  C  (In  l  )  ,  and 
orientation  angleD(in  rads)  of  a  cylinder  within  a  Couette  flow. 
D,(1  -20  1. 

Same  cylinders  and  liquids  properties  as  of  Fig.  5. 


Fig.  21.  The  vertical  velocity  A  (in  1/sec),  angular  velocity  B 


(in  rouds/sec)  ,  distance  from  the  wall  C  (ini),  and  orientation 


angle  D  (in  rods)  as  a  function  of  time.  Medium-air  of  p. 
-4 

»  1.718  X  10  poise;  particle  -  1  *  20  p»n.  •  b  *  0.5 
*2.23  g.  cm  channels  half  width  a  ■  0,5  cm; 

=0.15  ( 1 • sec)  '  . 
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'*•  Discussion 

The  results  shown  in  the  above  figures  should  be  discussed  while  keeping 
in  mind  the  behaviour  of  a  cylinder  under  creeping  motion  conditions  in  an 
unbounded  (far  from  walls)  fluid.  In  the  latter  case,  an  homogeneously  dense 
cylinder  which  settles  down  in  an  external  Held  devoid  of  orienting  forces 
moves  with  a  neutral  stability,  viz.  a  fixed  orientation  angle  (12).  In  our 
case,  however,  the  cylinder  is  definitely  subjected  to  a  torque. 

Vertical  Wall;  Still  Medium 

Here  one  can  observe: 

i.  A  continuous  change  of  the  orientation  angle  0  caused  by  the  wall- 
created  torque, 

ii.  A  continuous  development  of  a  drift  and  an  increase  of  the  displa¬ 
cement  D2‘ 

It  turns  out  that  the  sign  of  the  torque  depends  on  the  angle  ©  ;for  every 
initial  distance  of  the  particle  from  the  wall  there  exists  a  critical  angle  q* 
such  that 

Lf  <0  when  1 6  |  >  q  *  , 

L„  >  0  when  ie  I  <;  ©*  , 

and 

L  p  =  0  when  10  1= 

The  correspondence  between  theory  and  experiments  depicted  In  Figs.  5  through  7 
is  really  good,  which  produces  a  confidence  in  our  calculatlonal  method. 
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Horizontal  Wail;  Still  Medium 

Here  particles  motion  is  always  two-dimensional.  Its  characteristics  are: 

i.  A  continuous  decrease  of  the  sedimentation  velocity,  Vj,  as  the 
plane  is  approached, 

ii.  An  enhanced  wall  influence  for  small  , 

iii.  A  negligigble  effect  for  large  e0(  >  o.  “0  J 

iv.  An  increased  value  and  sensitivity  to  e0  of  the  angular  velocity 
as  the  particle-plane  gap  is  reduced, 

v.  A  constant  drift  velocity,  V2  ,  till  very  small  distances  from 
the  plane. 

It  is  interesting  to  compare  these  results  with  those  of  a  spherical  particle- 
horizontal  wall  situation(2)  where  a  similar  behaviour  is  observed. 

Inclined  Wall;  Still  Medium;  Two-dimensional  Situation 

From  the  results  for  the  inclined  wall  it  is  seen  that: 

\ 

i.  In  the  case  of  a  vertical  initial  orientation,  the  torque  is 
positive,  which  will  cause  the  cylinder  to  drift  from  right  to  the  left  (Fig. 12). 

ii.  In  the  case  of  an  horizontal  initial  orientation,  the  torque  is 
negative  but,  nevertheless, the  cylinder  will  drift  in  the  former  direction  too 
(Fig.  13). 

iii.  In  the  case  of  a  wal 1 -incl ination  of  a  little  less  than  75°  for  a 

vertical  initial  cylinder-orientation  and  of  a  little  less  than  60°  for  an 

horizontal  wall-inclination,  the  angular  velocity  ca  and  the  horizontal 

separation  02  respectively  reach  a  minimum  and  a  maximum  at  the  same  time. 

Thus,  the  particle  will  not  collide  with  the  surface,  which  is  of  a  great  slgnifi 

cance 


p 
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Here  one  should  note  the  aerosol  specificity  as  compared  to  the 

motion  within  our  liquid.  Since  both  the  mass  and  the  moment  of  Inertia  of 

2 

a  cylindrical  particle  are  proportional  to  b  ,  it  comes  out  that  the  thinner 
is  the  particle  the  greater  are  its  translational  (in  1/sec)  and  rotational 
velocities  attained  after  a  certain  time.  However,  if  one  returns  to  dimen¬ 
sional  quantities,  one  observes  that  the  overai 1  particle  drifts  in  aerosol 
and  1 iquid  systems, as  caused  by  a  plan  ar  wall,  will  still  be  the  same.  For  an 
aerosol  system,  the  influence  of  the  wall  is  manifested  already  at  distances 
of  about  50  1  (  2) . 

Motion  within  a  Flow;  Two-dimensional  Situation 

From  Fig.  20  it  can  be  seen  that  the  behaviour  of  our  cylinder  as  it  rotates 
in  the  gradient  flow  is  essentially  the  same  as  that  far  from  a  wall  ("Depo¬ 
sition  from  Flow",  Final  Report  DAERO-76-G-J4,  1978,  and  Gallfly  and  Eisner, 

J.  Colloid  Interface  Sci .  613,  320-37  (1979)).  The  wall  effect  is  shown  In  the 
decrease  of  the  second  maximum  of  the  vertical  velocity  V,  as  the  particle 
nears  the  horizontal  plane  (curve  A);  the  effect  is,  however,  over-shadowed  in 
the  angular  velocity  B  curve.  It  Is  noteworthy  that  the  (rotating)  particle  spends 
more  than  85$  of  its  time  in  an  horizontal  orientation. 

The  behaviour  in  this  case  Is  essentially  the  same  as  that  within  a  Couette 

flow. 


General  Conclusions 


i.  For  cylindrical  particles  there  is  a  pronounced  wall  effect  which 
may  be  of  a  great  importance  in  a  motion  within  small  cavities  like  those 
of  the  human  iung,  and  near  small  aerosol  collecting  obstacles. 

ii.  The  sedimentation  efficiency  of  the  particles  on  an  inclined  or 
horizontal  plane  is  decreased  due  to  the  effect. 

iii.  The  surface  concentration  of  the  part icles’  depos i t  on  an  Inclined  plane 
has  a  greater  dispersion  as  compared  to  that  for  which  the  wall  does  not 

exert  any  influence. 

iv.  The  collision  of  the  particles  with  a  vertical  wall  is  very  small  due 
to  their  "repulsion"  from  it;  only  particles  with  O  <  6  <^/2  will  have 

a  slight  chance  of  colliding. 

Finally,  extending  the  results  to  motion  near  curved  convex  surfaces,  one  can 
assume  that  the  smaller  is  the  angle  between  the  tangent  plane  and  the 
surface,  the  greater  is  the  (still  medium)  collision  efficiency  which  is 
defined  as  the  ratio  between  the  number  of  settled  particles  to  the  number  that 
would  have  sedimented  had  all  trajectories  been  vertical. 

V .  Applications 

The  slender  body  method  of  calculating  the  fluiddynamlc  force  and  torque 
on  the  particles  seems  to  be  more  useful  than  that  based  on  Oberbeck's  (1*) 
and  Jeffery's  (5)  equations  since  it  does  not  depend  on  any  limiting 
assumptions  concerning  the  flow  in  which  the  particle  is  Immersed.  For  example, 
any  numerical iy  (computer)  solved  field  of  flow  can  be  introduced  into  the 
slender  body  equations  where  It  is  necessary  tc  know  the  undisturbed  velocity 


(Eq.  Il9]  )  at  the  middle  of  each  of  the  N  divided  cylinder's  sections. 

This  has  been  performed  (with  no  wall  effect .however)  In  a  DFG  (Germany) 
supported  study  in  calculations  which  checked  quite  satisfactorily  with 
experiments  (13)- 


B.  BROWNIAN  DIFFUSION 


I .  Theoret i ca 1 
1 .  Method 
a.  General 


The  diffusion  of  aerosol  particles  In  air  is  generally  controlled  by  a 
Brownian,  a  turbulent  or  a  mixed  Mechanism;  near  surfaces,  where  turbulence 
subsides,  it  is  dictated  mainly  by  the  first. 

In  the  case  of  spherical  particles,  there  has  been  already  much  study  on 
the  numerous  aspects  of  Brownian  diffusion,  but,  very  little  has  been 
performed  with  regard  to  the  process  for  the  nonspherical  counterparts.  The 
difficulty  arises  from  the  general  non-isotropic  nature  of  teh  fluiddynamic 
interaction  between  particles  and  medium  which  makes  it  necessary  to  consider 
the  spatial  orientation  and  the  rotation  of  the  former. 

A  useful  way  of  treating  the  problem  Is  based  on  Brenner's  formal  structure  (14) 
in  which  one  assumes  a  six-dimensional  location-orientation  hyperspace  where 
the  general  number  concentration  Is  given  by  the  solution  of  the  conser¬ 
vation  -equation 
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Here  ,  the  translational  ,  J  ,  and  rotational,  T  ,  fluxes  -are 
respectively  expressed  (for  ax 1 “symmetrical  particles,  for  example)  by 


and 
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For  the  physical  applications  of  diffusion, where  one  needs  to  have  ensemble 
orientation-averaged  quantities,  It  is  useful  to  define  a  total  concentra- 
tion  c,  translational  flux  ITt  ,  averaged  translational  velocity  <u>  and  a 

-t 

diffusion  diadic  by 


■t 

The  total,  orientation-averaged  translational  diffusion  diadic  DT  is  the 
operational  parameter  to  be  plugged  into  the  (translational)  diffusion  equation 
of  the  aerosol  ensemble.  It  Is  dependent  on  the  orientation  distribution 
function  of  the  particles  which  Is  the  most  important  property  of  the 
ensemble. 
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In  Brenners'  (formal)  treatment  (14)  and  for  an  axl -symmetric  particles'  case, 
one  has  instead  of  [27J 


HC 
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where 
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C  u 


,and 


<r(r,e  ,*)  =  C  Cri)  f  Ofe,  +  ) 


[31] 


In  that  case,  one  obtains  for  the  orientation  distribution  density  f  the 
governing  equation 


+  Tt  '  *  X  (  r, 

which  takes,  for  example,  the  form 
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where  fc>  0  are  the  polar  angles  and  the  rotational  Peci&t  number  rPe 
is  defined  by 


V 

Pt 


[34] 


For  a  quiet  medium,  the  orientation  density  f  has  been  calculated  to  be 
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For  a  flowing  medium,  the  convective-rotation  equation  of  axisymmetric 
particles  for  example,  viz. 

|f  =  A„-f  -  r,  •  C*  «)  [3« 

has  been  (firstly)  solved  by  Peterlin  ( 1 6)  when  the  Couette  conditions  and  the 
pertinent  Jeffery  equations  for  &>  (5)  hold.  The  solution  shows  f  to  be  a  sole 
function  of  the  particle  Peclet  number  . 


b.  The  Aerosol  Specificity 

From  the  above  it  can  be  observed  that  the  (crucial)  density  function  f 

r 

for  axi-symmetric  particles  depends  on  Dx  and  V  .  The  latter  have  been 
calculated  In  the  case  of  a  1 iquid  suspension  by  classical  fluiddynamic 
methods  (17,5);  however,  in  the  case  of  aerosol  systems  this  Is  not  adequate. 

The  reason  for  the  inadequacy  stems  out  of  the  fact  that  in  the  very  Instance 
of  the  only  (small)  aerosol  particles  which  show  a  pronounced  Brownian  motion, 

the  medium  cannot  be  considered  as  continuous  any  more.  Thus,  both  the  rotational 
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diffusion  coefficient  and  the  gradient-rotational  velocity  of  the 

particles  have  to  be  calculated  from  first  principles  by  statist  leal -mechanics 
means. 

c.  Immediate  Aim 

It  was  intended  to  study  (for  the  first  time)  the  rotational  diffusion 
coefficient  of  cylindrical  or  discoidal  aerosol  particles  which  form  a  good 
geometrical  model  for  many  systems.  The  task,  has  been  accompl i shed, as  reported 
now  and  also  published  under  the  title  "On  the  Stochastic  Nature  of  the  Notion 
of  Nonspherical  Aerosol  Particles  III.  The  Rotational  Diffusion  Diadic  and 
Appl icat ions"by  Alfred  D.  Eisner  and  Isaiah  Gallily,  J.  Colloid  Interface 
Sci.  82.,  214-33  (1980. 

The  gradient  rotational  velocity  of  the  particles  In  a  simple  (though  not 
limiting)  shear  flow  has  been  calculated  too  and  presented  in  the  publication: 

"On  the  Stochastic  Motion  of  Nonspherical  Aerosol  Particles.  IV.  Gradient 
Convective  Rotational  Velocities  in  a  Simple  Shear  Flow"  by  Alfred  D.  Eisner 
and  Isaiah  Gallfly,  J.  Colloid  Interface  Sci.  1982  (In  press),  which  is  not 
included  here. 

Finally,  a  third  (summarizing)  article  is  in  preparation. 

d.  Medium,  Particle,  and  General  Governing  Equations 

Out  of  the  generalized  Einstein  equation 

»  *»/  |*r  (*  *  C  )  f 37] 
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between  the  translational,  rotational  or  coupling  particle  mobility  and  th« 
respective  diffusion  tensors,  and  from  the  relation  to  the  operating 
torque  u  exemplified  by 


Lf  =  -  H  .  «  ,  [38J 

one  notes  that  the  core  of  the  calculations  Is  imbedded  In  the  evaluation  of 
*  * 

The  (considered)  particles,  assumed  to  be  much  smaller  than  the  mean  free  peth 
1  of  the  gas  molecules  (  >>  i.  )  ,are  taken  to  translate 

with  a  linear  uniform  velocity  UM  with  respect  to  an  inertial  system 
of  coordinates  (Fig.  22),  and  rotate  with  a  uniform  velocity  u>  with 

respect  to  the  center  of  mass  M  In  a  quiet  medium  devoid  of  external  torques. 
Thus,  the  source  of  Lr  is  supposed  to  be  only  the  stochastic  bombardment 
of  the  gas  molecules. 

For  our  medium  (and  structureless  molecules),  the  distribution  of  the  molecu¬ 
lar  velocities  is  Maxwellian,  i.e. 

■dv?-  m  c*./,)*4  «.*?:->,<*.<  v“l.  m 

where  ( 1 ) ,  (2) , (3)  respectively  stand  for  X*  y*  ,1*  and  K  »  'wi/alrr. 
Concerning  the  Inertial  system  Ox  itself,  there  is  no  preferred  position  of 
Its  placement  in  space  since  the  latter  Is  assumed  (by  the  above)  to  be 
isotropic. 


Fig.  22.  The  external  system  of 
coordinates  0^  . 

The  particles'  surface  conditions  adopted  in  our  analysis  are  those 
of  Haxwell  (18)  according  to  whom  one  has  to  distinguish  between  a  diffusely 
reflected  fraction  t*  of  gas  molecules  with  a  Maxwellian  distribution 
conforming  to  the  surface  temperature  T,  and  a  specularly  reflected 
fraction  (  t  -  cc  )  for  which 


2  H  •  (  11  ■  Vt  ) 


I*»0] 
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From  the  point  of  view  of  an  observer  at  Oj,  the  velocity  of  a  gas 
molecule  relative  to  the  body-point  0^(FIg.  22)  is 

V,  =  Vt  -  Cjxf?  -  U  M  ,  [41] 


Thus,  the  number  of  these  molecules  of  velocity  between  vt  to  v .  *  )v, 
hitting  an  elemental  surface  ei  3  around  that  point  In  a  unit  of  time 
is 


rt("J  a -  -  N  (  W/rr)  ^[(y,  -  ox*  - 


•  exp[- k  iv*l 1  3  d3vt  ,  [42] 

the  momentum  transfer  at  that  time  is 

m  (v,-  «x  R-  UM)  -n'"1  d\  ol£  ,  C*3] 

and  the  torque  imparted  to  the  particles  in  the  fc  direction 
is 

(  Lf-  *)'  )  s  -  Nm  (h/rr)3//a  J{  [h.(*x  (v2  -  ox  R  -  u„))] 

£>  va 

X  o(5vx  dS  [44] 

which  Is  a  typical  governing  equation. 
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The  contribution  of  the  reflected  molecules  is  expressed  In  a  similar 
fashion  but  with  the  exchange  of  h  and  N  of  the  velocity  distribution 
of  the  incoming  molecules  with  h-r  and  N  of  the  diffusely  reflected 
ones  as  determined  by  the  requirement  that  (for  a  steady  situation): 


i.  The  net  flux  of  molecules  at  each  particles  surface  element  is  zero; 

ii.  The  net  total  energy  transport  at  that  element  Is  zero,  too. 


The  total  torque  imparted  to  the  rotating-translating  particle  can  be  put  as 


l,  ■  w;w>  = 


.  cl,.  »)”""]  .  [«] 


MCcO- 


2.  Calculational  Procedure 
a.  A  Spherical  Particle 

This  has  been  chosen  as  a  test  case  for  the  proposed  method.  Taking  here 
the  "natural"  coordinates  system  »  "i#  and  (Fig.  23)  •  one  has 


R  =  t~  (  Sim  &  cos  $  i  Sin®  si-rfi  +■  CO» 6  k)  , 


[46] 


and 


n  =  siwacestl  ■¥  s'i>i  &  srri  Q  j  +  cos  6  k 


[47] 


[48] 


1 
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where  A,  the  transformation  matrix  is 


Fig.  23.  The  spherical  coordinate  system 


Due  to  the  space  isotropy,  it  is  possible  to  write 


cj  s  wk 


and 


(  si n  S'  i  +  co*  M'  h  ) 


A  sample  procedural  way  for  calculating  (  .  |?  ) 

is  given  here  as  follows: 


which 


Also  , 


.  The  integrand 

I,  »  k  •  [  R  X  (V2  -  c JkR-UJ]  [(vrU*S  -  U*)*™] 

=  r  «in6  vaV  _  v/  -  r-^-co  im*e  v/ 

+  y 1  cO  UM  COS  0  ilyt2S  . 

.  Integration  boundaries 

,  r 

4o  arid  so  }  as  rise 

v  «  *  v.r  =  cos  e 

gives 


-  o*  4  v*.  4  UM  cos  0 

-  <*>  4.V4  t  <*> 


i  i  i .  Final  equation 


(L  •k/"’  ■=  -ti-m  (K/rr)^  /  f  f  [  f  I, 

0  S  .» 


/  e'  ^  Yi9  +  v/1)  |  J |  ot  v»r  alv»* <J[v*  6  «!♦  [56} 


T” 
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and 


where  B,  the  transformation  matrix,  is 


Placing  the  02(external)  coordinate  system  (Fig.  24)  so  that  the  polar 
axis  of  the  particle  will  coincide  with  I?  ,  and  assuming  an  arbitrarily 
directed  angular  velocity  c»>  =  co,  i  +  +  co, 

one  can  carry  out  a  piecewise  calculation  of  the  total  torque  Imparted  to 
the  envelope  and  each  base  of  the  particle.  By  order  of  magnitude  conside¬ 
rations  (APPENDIX  2),  it  can  be  shown  that  the  ordinary  translational 
velocity  of  an  aerosol  particle,  \JM  ,  is  negligible  with  respect  to 
and  cj  x  f?  { Eq .  [44]  ).  Consequently,  it  becomes  possible  to 

express  the  integrand  of  [44J  for  the  molecules  which  hit  the 


particles  envelope,  for  example,  as 
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I  =  i  •  [  f?  X  (Vj  -  <JX  )?;j  T(vi  -  at  x  R).  nj  =  -  i  *>**  V 

+  (-  *''"A4>  -oj,Zi<iV«a^+  <w*  tac©K*>  s«*  *  +  wa*aCo«*  Sm* 

+.  c*)3  j ? cos  4>)  <4*  +  (  cj,  $ z  -  Wj ; Z  st«  ♦  «*  ♦)  4  +  C-  «*>,  *  **>♦  cm * 

+  Mi  i2C0Sa  $)  Vj.^  +  ?  si-M  4>  V2*  Va*-  *  V»T  v/  -  W,H3S !•»»♦ 

+  tt,1?1  Z- cos2**.*  #  -  co2V*  to*4**?**  -  ^Vz*'**  +  2  W.WxJ^t 
co6  *  si-nT4>  +  CJ,  coa  i!J  cos  $  +  «,  w3  5  2 2  00s  4>  s  i'-n  4>  -  ca*  co3  j  Z4 

cos1*  .  [63J 


The  integration  limits  in  the  (governing  equation)  [44]  are  set  to  be 
V,  >  n  =.  vj  +  co,  ;  iw  ♦  -  cc>i  Z  cos  4>  —  0  f  which  leads  to 

-00  £;  v/  £  cj2  2  004  4>  -  «0,  2  n  -n  ♦  , 

4>  ^ 

and  -®  -  V,  j  Vi  4  >» 


Final lly ,  in  this  case,  the  torque  can  be  written  in  the  form 
(lf.  i)c0,°  =  -^  (h/TT)V*;/  f  f  I, 

F  c  Sf6  V*  Vf  xt 

d  v/  o(  ^  d  Vi*  d  £  [64] 

by  the  above,  and  dJ>=?^ao*^ 


where  the  velocity  limits  are  given 
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For  the  case  of  ^  Cjiz<< 1  which  almost  always  occurs  (APPENDIX  2),  and 
when  utilizing  certain  integral  values  (APPENDIX  2),  one  obtains  in  the 
typical  procedure  shown  here  that 


(i*) 

(u  F  •  t  )c  =  -  N  1»?  (  w/rr)  l [- 

$  «  h  J 


I  /  i.  .2 . 

-JJJ  (-  «.  J  *•*»  0 


CJ,  4acos1<fi  +  coj.  f2  OOJ  #  si'n  ♦  -  coj*3c«*a4  tiw  ♦  +  cjjS?  Co«f) 


+  (co.Va  cos14)i^0+ 

+  CO,  COi  Z3C0S,<J>  +  OO,  Cd^  ?a2Si'<y»*4>  ■+  CJ,  COa  S^Z  Cos#  Si'Wi#  1  CO ,&*,  jr1 

0>s#  Si'-n  f  —  coacyj  co^H  cos#-  aj,i*-s»n  *)+£>(-)],#  [651 

where  0(...)  denotes  higher  powers  of  co;  c o^'  which  can  be 

neglected  (APPENDIX  2). 

I  I .  Resul ts 


1.  The  Spherical  Particle 

The  total  torque  on  the  spherical  particle  for  to  *  co  k  and 
is  found  to  be 


(lf  •  K)(t#+)= -o((  (■ 


_4_ 

3 


IS 


h  U*  ) 


[66J 
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X  C-|-  -  -i|~  ( ' +  ^  u„a  ]  [fi93 

which  becomes  identical  with  [66J  as  =  O  . 

2.  The  Cylindrical  or  Discoidal  Particle 
The  total  torque  on  these  particles  Is 
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0-f'  l)c  r c ~zf- ■*■  ) 


f )  *  «  [  (4—  -fr>i '("•**  “i> 


8  CJ, 


r  -il  i* 

^  16  +  80 


8  L 


+  -hr)'!  +  OC--0  t70] 


which  can  be  written  to  a  very  good  approximation  (APPENDIX  2)  as 


+  - ;  P-  i?+  ^P8)! 


[71] 


where  f*  -  §  /  *- 


Similarly  to  Eq.  71  , i t  is  possible  to  get  an  equation  for  the  torques 

imparted  to  the  rotating  aerosol  particle  along  its  other  axes;  these  come 
out  to  be 


( L  r  •  i  )  rv  -  CJi 


[72] 


and 


(Wr  -W)  ~  -  co3 


C73] 


Thus,  the  (vector)  torque  can  be  written  as 


where  -fL  =  M 


L  h  =  -  il.  •  cj 


is  a  symmetric  diadic  whose  eigenvectors  lie 


along  the  symmetry  axes  of  the  particle.  A  similar  deduction  can  be 
reached  from  Landau's  and  Lifshltz's  (19)  reasoning  based  on  Irrever¬ 
sible  thermodynamics  principles. 


a.  The  rotational  diffusion  coefficient 


From  Eq.  [38]  one  has  that 


(LP.  i)  ^ 


which  gives,  if  one  uses  the  kinetic  theory  relationship  Uy*x{  jx  (n/T) 


rD  =  _ WT  tt-n _ 

"  11  'f4  l3  l  (  1  *  T?  <  A?*)  -  Ap. 

where  Kn  =  L/f  . 


Eq.  [76]  should  be  compared  to  that  of  Edwards  (20) 


8  CL'/3)(  -Jy-j-) 


[77] 
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as  is  presented  in  Fig.  25. 


Fig.  25.  The  rotational  diffusion  coefficient  of  Cylincrical 
particles  in  air.  Af,Bj,Cj  -  values  of  this  study; 
Aj.B,C  -  Edwards'  values. 


111.  Discussion 

It  must  be  re-noted  here  that  our  calculations  apply  to  aerosol 
particles  whose  dimension,  in  any  direction,  is  much  smaller  than  the 
mean-free  path  1  of  the  gas  molecules  (*0.069  at  normal  conditions), 

and  to  a  quiet  medium.  For  larger  particles  within  a  flow,  one  has  to 
recur  to  the  Boltzmann's  Integro-dlfferentlal  equation  of  the  velocity 
distribution  whose  (closed  form)  solution  is  not  known  as  yet. 
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The  results  of  0A  themselves  show  that: 

i.  The  values  of  the  coefficient  are  of  the  order  of  10^-10^  sec 
whicli  makes  the  typical  relaxation  times  for  Brownian  rotation  phenomenta, 

-  =  '  i_  negligibly  small. 

ii.  The  calculated  values  are  greater  by  a  factor  of  about  30 

than  those  obtained  according  to  Edwards  in  classical  fluiddynamic  means,  which 
is  reflected  by  the  rotational  Peclet  number  assigned  to  a  diffusional  con¬ 
dition  of  orientable  particles  and,  hence,  by  the  orientational  distribution 
densi ty. 

IV.  Appl i cat  ions 

As  an  example  for  the  use  of  in  a  real  problem,  it  is  possible 

to  take  the  diffusional  deposition  of  cylindrical  aerosol  particles  on  an 
infinite  plannar  surface  and  in  a  quiet  mediuw. 


3  .  General 

Employing  the  cartesian  coordinate  system  shown  in  Fig.  26,  one  can 
write  down  Eq.  [27] in  the  form 


§f  =  [  o,  *  (X-Vx)  cos4*!  +  Z\*  (*!>»-' 


Si7,as  s;«3*J  -ll§_  +  c  +  -Vl)  “***1 

v  3 


IV 


4  2  (d„  -  )[sir?*6  sin  *  CO*  *  — ^ — £ -  +  si 

"  ax,  3 


r — £ -  +  S  in  0  cos  d  Sin* 


•SC 


s<n0nos0  nos* 
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Fig.  26.  The  cartesian  system  of  coordinates  used  in  the 
diffusion  analysis. 

where  (|  }  *33  x  are  the  translational  diffusion  coefficients  of  the 

particle  parallel  and  perpendicular  to  its  polar  axis,  respectively.  Now, 
if  it  is  assumed  that  particles'  orientation  density  is  Independent  of 
thei r  location  viz . 


<r(r,#,t)=  c  (*.  ,  Xj,  x,)  ♦  .  -O  ,  f79l 

then,  instead  of  [783  ,  one  obtains  the  two  equations: 


34 /M  =  *>±1 


[80) 


whose  solution  is  given  In  Eq.  [35>],and 
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-i&-  =  z  K  +  C  x>H-Vj  •«a«co»a4ir£4+c  i>„-^) 

o  ir  1 


a1/*  ♦  ,  4  * 

,  r  ^  ^  f  n>  _  -r 


Si^  <mV]  -|X  ♦  C  !dx  +  (  *>„  J  Mi**  ]  -fe 

O  Xj  3 


2  ^  P„  -  J3a)  [  »ir2e  sin  4>  coi^^)  — — — - —  +  sm*  cos*  sir)  ^ 

ax, 


«fs,+  s‘re  “‘s  ^*^nx ,3  ■ 


At  this  stage,  the  orientation  averaging  being  defined  by 


*n  tt 

,fj  3  l  ( 6>  ttfsine  de  d  <P 

rf^  tir< a  ole  d$ 


one  reaches  the  final  equation 


N  3t  ' 


J  S  "It*  ^ 

3[  j"  J-  (e^t  ,4)  S'”G  Jed<P 


=  [  x>x  +  Cv„  -+X)x )  <  sin1©  cosJ<*>>]  -|-~r  +  C  C-t>„- x>x) 


<SiVe  si-na4»>J  1—  +[ii  +  (  J5A+(2>H-^x)<COSae>]-^-  [&2] 


since  the  mixed  terms  9*C  /  SX;  3  Xj  drop  out. 


b.  Specific 


In  the  specific,  above  diffusion  case  where  it  is  taken  that  concen* 
tration  gradients  exist  only  along  #» 


If  - 


t 
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where  the  orientation-averaged  translational  diffusion  coefficient  ob  (t) 


i  s 


cZ>  Ct )  =  +  c  »„  -\)<siy,ze  .  [84] 

Peforming  the  orientation-averaging  of  the  function  si-*1®  nV  ♦  one 
gets  that 


6bCt)=X  +  (x>ll-K){-f"C»-^c-6*x+)]}.  tsa 

Concerning  the  solution  of  [83]  ,  it  is  advantageous  to  define  a  (new) 

vriable  P  by  (21) 

t 

P  =  /  oZ>Ct’)d  +  '  [86] 

0 

whereby  one  obtains  that  (for  c  =  c0  at  t  »  o  and  c  *o  a.t 

xa  =  o  and  t  >  o  ) 

c  =  Co  $  C  )  *  c*  *  C  2  ( 3  »  t87j 


§  being  the  Error  Function. 

The  total  number  of  deposited  particles  M^after  time  t  will  be  then  (22) 


-  >/  - 


m+  =  ict  - : - ) 

■  v  rr  ' 


c .  Results 


The  number  M  -  vs.  time  is  presented  in  Fig.  27  for  the  two 


A.  When  *D  and  +  DX  are  calculated  according  to  classical  continuum 
theory,  viz. 


X  it  ft.  L  [  '/< f-n  (i/p's  - 


^fTT  f*  l  [•/<  1*  C'/P)  +  o.r>3 


and  Edwards'  value  of  X>  A  (  5  ^j.  )  's 


assumed. 


B.  When  these  coefficients  are  given  by  Dahnekes  molecular  regime  equatlons,vlz. 


D  „  =  — ' 


_ bJL±J2 - . - [91] 

tt  f  J  [  q  +  of  <” ;  +  (i/P)-  4-rr?A>3 


- SL2_l23 -  [92] 

^  ^  S  C (1/1»)  +  «< ' +  f '/P)  -  C  t 


and  our  value  of  X>  ,  Is  taken. 


58 


Also,  ^ _  is  calcualted  for  case  C  of  a  volumetrical ly  equivalent 
spherical  particle. 

From  Fig.  27  it  is  seen  that  the  values  of  case  B  are  significantly  higher 
than  those  of  case  A  but  that  those  of  case  C  are  higher  still,  which  is 
resonable  if  one  remembers  that  diffusional  particle  mobility  is  a  property  of  the 
surface  of  tne  latter  and  that  the  surface  of  the  volumetrical ly  equivalent 
sphere  is  the  smallest  of  all  forms. 


it  should  be  noted  that  tn  the  st i 1 1 -medium  diffusion  discussed  here 
is  essentially  independent  of  the  (very  large)  value  of  o£?  (t)  for  all 
tal  times  of  interest.  The  value  of  D  t  will  matter  in  difusion 
where  the  rotational  Pec  let  number 
count . 


Pep  (  =  leal  /  x  j 


ex per imen- 
in  flow 
does 


a»J 

012 


l 


Fig.  27.  The  total  number  vs.  t. 
A,  case  A-,B,  case  B  ;  C  ,  case  C  cited  above  . 


f  =  O.  0  015  ,  L  ~  0.1  T=  fi-  1»q.  22  *  < 0'{  poise 

1=  6.  6  9  x  IP-*  cm. 


5? 


C.  APPENDIX  1,2 

1 .  Append i x  1 

a .  Green  functions 

General lv,  expressing  the  disturbance  due  to  the  slender  body 
by  (Fig.  1) 


where 
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fj-CrO  -  $  J*  , 

Cr')  =  ^  f  j  <  r‘>  ?  ^  ^  +  ** C  ' 

s  i  nee 

*->•  =  r1  ♦  ffr',  O  , 

b.  Definitions 

£  s  C  ,’  c*sfl  +  Ca  i'y' e  j  C2  =  -C'tinS  *  C'2  c*s  0  ,  Cs=  C3  # 

M((  =  2  tos^fl  +  S  iVd  Ml2  =  M2I  =  -  si-ne  CM0J  M„  =  MJ(  =  0; 

Mns  2  s.-n-1©  •*  CosJ<9  M23  =  M  a*  =  O  }  M  M  =  1  . 

TT,,  ss  K  coi*0  +  K22  siVS  j  TTJJt  =  TTj,  =:  K)SiV,e  «)*  ^  ^ 

TT„  =  !<2j  S"T7  0  3  rrJ3  =  S .••*»*«  +  K»  Oos3©^  TT-25  =  K2,C03«, 

tt81  =  Kw  sw  e  ,  ttJ4  =  kt32  c-oi  &  ,  TT,,  =  i<„  . 

r„  - ■■  ( f?  /  ?*)  si-n'e  ,  iu=  it,  =  (tl/f1)  cose, 

I,3s  I31  =  12  = 

i  '  =  sr*a  ,  i,2  «  =  sf-n  e  cos  e  ,  x  (3'  =  o  ,  x«  *  001 

i«  =  xi  =  xa:  =  o. 

X(t"  =  00**0  ,  X,2"=  Ia,"  =  -si*0  cot^ 2  Xa2  =-6‘^e,  1i»"  I»i  ~ 
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and  writes  £  f  7J  and  [18]  in  terms  of  these  transformat  ions . 

Thus,  u^w(*jx’)  - »Aij  (v'^  and  Eq.  [2CQ  is  expressed  as 

u'V  =E  r.r>  ir 

2T* 

where  -  N  (  i  -1 )  +  '  =  hJ^-O  +  tr'  , 

=  -y-  m  c  jf, r’ ) +  ° r  r * 2r ■  >  'r(r > jt')  =  N CjTjr’^+>v^y’'>  **r 3"a r- 


2.  Appendix  2 

\  in 

a.  Order  of  magnitude  of  h  Um 

Taking  a  typical  case  of  =  0  (i)  and  calculating  h  to  be 

5-^7  *  ,o'  10  sec  2  cttT*  for  T=  2<*8.ifc*K  and  a  i  r  molecules , 

one  obtains  that 

h'/lUp1  =  0  (  lO"S)  ,  which  is  much  smaller 

than  unity. 

b.  Order  of  magnitude  of  V> r.  cot 

Taking  L  =  0.1  x  I0'*1  tw  ,  <j=o.0O75'yl0  ^  c-m.^  =  1^.  Ctrt  s 

and  air  molecules  at  T  *  118  .  U  '  K  ,  t*  *  1  8  *  ■ 1 2  *  »(T*  po  is<,  ^  cw 


and  estimating  to,  by  the  law  of  equiparti tion  of  energy,  viz. 


t 

T 


It  <*>«*>  = 


te  t 


fC.3] 


one  gets  that 
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D.  NOMENCLATURE 

a  -  constant  of  systems  geometry;  half-width  of  flow  channel 
A, 8  -  transformation  matrices  defined  in  text 

b-  characteristic  equivalent  cross-se ' t i on  radius  of  (slender)  particle 

b  -  local  equivalent  cross-section  radius 

S;  ,  B ,  B  -  constants 

c  -  particles'  total  number  concentration,  c  -  initial  value 

C  -  vector,  defined  in  text 
D  -  distance,  defined  in  Fig.  2 

t  e 

0  ,  D  .  D  -  t rans 1  a t i ona 1 ,  rotational  and  coupling  diffusion  diadics, 
respect iveiy 

Oj.  -  rotational  diffusion  coefficient  around  a  mid-diameter 

V  i  V  ~  translational  diffusion  diadic  components  parallel  and  perpendicular  to 
polar  axis,  respectively 
0,  ,Dj  -  distances,  defined  in  Fig.l* 

Jb  -  4  Cb  -  tux)  <  s  v<9  r,  > 

e  -  unit  vector  along  symmetry  axis  of  particle 

f  -  orientation  distribution  density 

-  Stokeslets'  distribution  density;  f^te(v)  -doublets'  distribution 

dens i ty  etc . 

-j-  (r"  -  drag  force  per  unit  area  of  particle 

FP  -  external  force;  Ff.  -  fluiddynamic  force 

g  -  acceleration  of  gravity 

h  -  particle-wall  distance;  also,  m/2kT 

If 

I  ,  It  •  moment  of  inertia  of  particle 
I  ,  I',  I"  -  diadics,  defined  in  text 
I,  -  intergrand,  defined  in  text 

j  -  rotational  flux  of  ax i -symmetr i c  particles 

J  -  particles'  flux 

I J I  -  transformation  Jacobian 

k  -  Boltzmanns  constant 

k,  ,  k  -  constants 

Kvi  -  Knudsen's  number 

1  -  semi-length  of  slender  particle 

1  -  mean  free-  path  of  gas  molecules 

L  -  length  of  cylindrical  particle 

U  e  -  external  torque; 


fluiddynamic  torque 
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m  -  mass  of  particle 

M  -  total  number  of  deposited  particles  in  time  t 

M  -  diadic,  defined  in  text 

-  mobi li ty  diadic 

Tt  -  unit  vector  normal  to  dS 

n,  n'  -  running  indices 

~  '  ,  v  -  number  of  molecules  hitting  unit  surface  per  unit  time 

‘  -  density  of  Maxwell’s  velocity  distribution 

N  -  average  number  of  gas  molecules 

P  -  defined  in  text 

r. ..  -  rotational  Peclet  number 

'Pef  -  particle  rotational  Peclet  number 

-  associated  Legendre  function 

K  -  typical  length  of  particle 

v)r-'|i r".  r*  -  distances  defined  in  Fig.  1  and  2 

r  -  location  radius  vector  ,  v  -  rnd'us  o4  sphere 

^  -  radius  vector  from  geometric  center  of  particle  to  a  point  on  its 

surface. 

S  -  particles  surface;  ot S  -  surface  element 

t  -time 


T  -  Temperature 

cl,  (r  r',  -  etc.  -  Green  functions  associated  with 
UM,a  -  translational  velocity  of  particle 
v  Cr)  -  disturbance  velocity  of  slender  body 

vX;  v,  -  velocities  of  gas  molecules  in  the  0X  system  and  relative  to  dS 
respective1y;v’  -  velocity  of  refected  molecules. 

V  -  translational  velocity  of  slender  body 

v  Cr)  -  double-layer  potential 

X  t  X1  -  position  distances  on  body's  axis  (Fig. 1 ) 

X,  ,  Xj  ,  X5  -  coordinates 
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Greek  letter-' 


-*  -  friction  of  gas  molecules  reflected  diffusely 

s  -  ? 

-  running  indices 
t  -  Kronickers1  delta 

A  -  1  /  N 

-  v-  !  evi-Civita  pc-rr.iul  at  ion  tensor 
'  -  fluid's  ror t i ci ty 

*  -  angle,  defined  in  text 

<■.*  -  polar  angles 

-  vector,  defined  in  text 
■  -  dynamic  viscosity  of  gas 
TT  -  diadic,  defined  in  text 
f  -  position  vector,  defined  in  Fig.  1 
?  -  gas  density;  cy,  nr.  dees  raou'us 
°c  -  particles'  density 

a  -  number  concentration  in  location-orientation  hyper-space  • 

a  1  sc'  j  i_  (.y, ,  2  ..  o  J 

4>  -  o-ientation  (pseudo) vector 

$  -  Error  Function 

N’  -  angle  between  \Jt.  and  co 

ca-  rotational  velocity  of  particles;  co,  ,  co2  ,  co  3  -components 
A2.  =  m  ' 

Subscripts 
M  -  center  of  mass 
o  -  initial  value 
c  -  cylinder's  envelope 
Supe r  s  C r i o t  s 

(in)  -  hitting  molecules 

c*  -  t,r,c 

(r)  -  reflected  molecules 

(tot I  -  tCtrl 
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